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PARTIALLY ORDERED SETS IN MACAULAY2 

DAVID COOK II, SONJA MAPES, AND GWYNETH WHIELDON 

Abstract. We introduce the package Posets for Macaulay2. This package provides a 
data structure and the necessary methods for working with partially ordered sets, also 
called posets. In particular, the package implements methods to enumerate many commonly 
studied classes of posets, perform operations on posets, and calculate various invariants 
OQ \ associated to posets. 
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q ; Introduction. 

A partial order is a binary relation ■< over a set P that is antisymmetric, reflexive, and 
transitive. A set P together with a partial order ^ is called a poset, or partially ordered set. 

Posets are combinatorial structures that are used in modern mathematical research, par- 
ticularly in algebra. We introduce the package Posets for Macaulay2 via three distinct posets 
c^ ' or related ideals which arise naturally in combinatorial algebra. 

We first describe two posets that are generated from algebraic objects. The intersection 
semilattice associated to a hyperplane arrangement can be used to compute the number of 
unbounded and bounded real regions cut out by a hyperplane arrangement, as well as the 
dimensions of the homologies of the complex complement of a hyperplane arrangement. 

Given a monomial ideal, the lcm-lattice of its minimal generators gives information on 
the structure of the free resolution of the original ideal. Specifically, two monomial ideals 
with isomorphic lcm-lattices have the "same" (up to relabeling) minimal free resolution, and 
the lcm-lattice can be used to compute, among other things, the multigraded Betti numbers 
0i,h(R/M) = dim^ Toy i t b(R/M,k) of the monomial ideal. 

In contrast to the first two examples (associating a poset to an algebraic object), we then 
describe an ideal that is generated from a poset. In particular, the Hibi ideal of a finite 
poset is a squarefree monomial ideal which has many nice algebraic properties that can be 
described in terms of combinatorial properties of the poset. In particular, the resolution and 
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Betti numbers, the multiplicity, the projective dimension, and the Alexander dual are all 
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nicely described in terms of data about the poset itself. 

Intersection (semi)lattices. 

A hyperplane arrangement A is a finite collection of affine hyperplanes in some vector 
space V. The dimension of a hyperplane arrangement is defined by dim(^4) = dim(V), and 
the rank of a hyperplane arrangement rank(^4) is the dimension of the span in V of the set 
of normals to the hyperplanes in A. 

The intersection semilattice C(A) of A is the set of the nonempty intersections of subsets 
of hyperplanes f] n&A / H for 7i G A' C A, ordered by reverse inclusion. We include the empty 
intersection corresponding to A' = 0, which is the minimal element in the intersection meet 
semilattice G C(A). If the intersection of all hyperplanes in A is nonempty, f] nGA H ^ 0, 
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then the intersection meet semilattice C(A) is actually a lattice. Arrangements with this 
property are called central arrangements. 

Consider the non-central hyperplane arrangement A = {Hi = V(x + y),H.2 = V(x),H-3 — 
V(x — y), %4 = V(y + 1)}, where % = V(£i(x, y)) C M 2 denotes the hyperplane H t of zeros 
of the linear form £i(x,y); see Figure l(i). We can construct C{A) in Macaulay2 as follows. 

il : needsPackage "Posets"; 

i2 : R = RR[x,y] ; 

i3 : A = {x + y, x, x - y, y + 1}; 

i4 : LA = intersectionLattice(A, R) ; 

Further, using the method texPoset we can generate EKTgX to display the Hasse diagram of 
C(A), as in Figure 1(h). 





(i)A 




(ii) C(A) 



FIGURE 1. The non-central hyperplane arrangement 

A = {Ux = V(x + y),H 2 = V(x),H 3 = V(x - y),H 4 = V(y + 1)} 
and its intersection semilattice C(A) 

A theorem of Zaslavsky [Za] provides information about the topology of the complement 
of hyperplane arrangements in M n . Let fi denote the Mobius function of the intersection 
semilattice C(A). Then the number of regions that A divides IR n into is 

r{A)= Yl IM0,x)|. 
xeC(A) 

Moreover, the number of these regions that are bounded is 

b(A) = HC(A)U 1)1 

where C(A) U 1 is the intersection semilattice adjoined with a maximal element. 

We verify these results for the non-central hyperplane arrangement A using Macaulay2: 

i5 : realRegions(A, R) 

05 = 10 

i6 : boundedRegions(A, R) 

06 = 2 

Moreover, in the case of hyperplane arrangements in C n , using a theorem of Orlik and 
Solomon [OS], we can recover the Betti numbers (dimensions of homologies) of the comple- 
ment A^yt = C n — VJA of the hyperplane arrangement using purely combinatorial data of the 
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intersection semilattice. In particular, A4a has torsion-free integral cohomology with Betti 
numbers given by 



Pi(M A ) = dim c (^(M A )) = J2 ^ 5 > 



xGC(A) 
dim (x)=n— i 

where //(•) again represents the Mobius function. See [Wa] for details and generalizations of 
this formula. 

Posets will compute the ranks of elements in a poset, where the ranks in the intersection 
lattice LA are determined by the codimension of elements. Combining the outputs of our 
rank function with the Mobius function allows us to calculate 0o(Ma) — 1> Pi{M-a) — 4, 
and /3 2 (Ma) = 5. 
i7 : RLA = rank LA 

o7 = {{ideal 0}, {ideal (x+y), ideal (x), ideal (x-y), ideal (y+1)}, 
{ideal(y,x), ideal(y+l,x-l) , ideal (y+1 ,x) , ideal (y+1, x+1)}} 
i8 : MF = moebiusFunction LA; 

i9 : apply(RLA, r -> sum(r, x -> abs MF#(ideal 0_R, x))) 
o9 = {1, 4, 5} 

LCM-LATTICES. 

Let R = K[xi, . . . ,Xt] be the polynomial ring in t variables over the field K, where the 
degree of Xi is the standard basis vector e^ G Z*. Let M = (mi, . . . ,m n ) be a monomial 
ideal in R, then we define the Icm-lattice of M, denoted Lm, as the set of all least common 
multiples of subsets of the generators of M partially ordered by divisibility. It is easy to see 
that Lm will always be a finite atomic lattice. While lcm-lattices are nicely structured, they 
can be difficult to compute by hand especially for large examples or for ideals where Lm is 
not ranked. 

Consider the ideal M = (a 3 b 2 c, a 3 b 2 d, a 2 cd, abc 2 d,b 2 c 2 d) in R = k[a,b,c,d]. Then we 
can construct Lm in Macaulay2 as follows. See Figure 2 for the Hasse diagram of L M , as 
generated by the texPoset method. 

ilO : R = QQ[a,b,c,d] ; 

ill : M = ideal (a~3*b~2*c, a~3*b~2*d, a~2*c*d, a*b*c~2*d, b"2*c~2*d) ; 

il2 : LM = lcmLattice M; 

Lcm-lattices, which were introduced by Gasharov, Peeva, and Welker [GPW], have become 
an important tool used in studying free resolutions of monomial ideals. There have been 
a number of results that use the lcm-lattice to give constructive methods for finding free 
resolutions for monomial ideals, for some examples see , [CI], [PV], and [Ve]. 

In particular, Gasharov, Peeva, and Welker [GPW] provided a key connection between 
the lcm-lattice of a monomial ideal M of R and its minimal free resolution, namely, one can 
compute the (multigraded) Betti numbers of R/M using the lcm-lattice. Let A(P) denote 
the order complex of the poset P, then for i > 1 we have 

f3 hb (R/M) = dimH t _ 2 (A(6,b);k), 

for all b E L M , and so 

Pi{R/M) = J2 dimH^ 2 (A(6,b);k). 

beL M 
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x a 3 b 2 c 2 d 




a 3 b 2 cd 



a 3 b 2 d 



Figure 2. The lcm-lattice for M = (a 3 b 2 c, a 3 b 2 d, a 2 cd, abc 2 d, b 2 c 2 d) 

These computations can all be done using Posets together with the package Simplicial- 
Complexes, by S. Popescu, G. Smith, and M. Stillman. In particular, we can show that 
fii,a 2 b 2 c 2 d — for all i with the following calculation. 

il3 : Dl = orderComplex(openInterval(LM, 1_R, a~2*b~2*c~2*d)) ; 
il4 : prune HH(D1) 



ol4 = -1 

1 









ol4 : GradedModule 

Similarly, we can show that Pi^tfcd — 2. 
il5 : D2 = orderComplex(openInterval(L, 1_R, a~3*b~2*c*d)) ; 
il6 : prune HH(D2) 
ol6 =-1:0 

2 
: QQ 
ol6 : GradedModule 

HlBI IDEALS. 

Let P = {pi, . . . ,p n } be a finite poset with partial order ^, and let K be a field. The 
Hibi ideal, introduced by Herzog and Hibi [HH], of P over K is the squarefree ideal Hp in 
R — K[xi, . . . , x n , 7/i, . . . , y n ] generated by the monomials 

ui =W Xi Y[ ^ 

PiGl Pi0 

where I is an order ideal of P, i.e., for every i G / and p G P, if p ^ i, then p e /. Nota 

bene: The Hibi ideal is the ideal of the monomial generators of the Hibi ring, a toric ring 

first described by Hibi [Hi]. 

117 
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P = divisorPoset 12; 
HP = hibildeal P; 
HP * 



ol9 ={xxxxxx , xxxxxy , xxxxyy , xxxxyy , xxxyyy , 
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012345 012345 012435 012345 013245 
xxxyyy , xxyyyy , xxyyyy , xyyyyy , yyyyyy} 
012345 021345 012345 012345 012345 

Herzog and Hibi [HH] proved that every power of Hp has a linear resolution, and the 
i th Betti number (3i(R/Hp) is the number of intervals of the distributive lattice C(P) of 
P isomorphic to the rank i boolean lattice. Using Exercise 3.47 in Stanley's book [! ], we 
can recover this by looking instead at the number of elements of C(P) that cover exactly i 
elements. 

i20 : betti res HP 

12 3 
o20 = total: 1 10 12 3 




i21 
i22 
i23 
i24 
o24 
i25 
o25 



1 . . . 

. 10 12 3 
LP = distributiveLattice P; 

cvrs = partition(last, coveringRelations LP); 
iCvrs = tally apply (keys cvrs, i -> #cvrs#i) ; 
gk = prepend(l, apply (sort keys iCvrs, k -> iCvrs#k)) 
{1, 6, 3} 

apply(#gk, i -> sum(i . .<#gk, j -> binomial (j , i) * gk_j)) 
{10, 12, 3} 



Moreover, Herzog and Hibi [HH] proved that the projective dimension of Hp is the Dil- 
worth number of Hp, i.e., the maximum length of an antichain of Hp. 

i26 : pdim module HP == dilworthNumber P 
o26 = true 
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